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abstract 


The.  problem  of  approximating  the  terms  of  a  sequence  was 
investigated  by  M.  G.  Darboux  in  1878.  Since  then  the  theory 
of  asymptotic  expansions  has  been  greatly  extended.  The 
purpose  of  this  thesis  is  to  obtain  the  Darboux  result  by  more 
modern  means  and  in  a  form  more  consistent  with  present  theory. 
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INTRODUCTION 


The  concept  of  asymptotic  expansions  originated  with 
H.  Poincare"  in  1886*  The  usefulness  of  these  expansions  has 
resulted  in  numerous  generalizations  of  this  initial  definition. 
One  of  the  most  recent  and  elegant  of  these  definitions  is  due 
to  A.  Erdelyi  and  is  discussed  in  chapter  one,  A  specialization 
of  this  definition  is  adopted  for  the  purpose  of  the  thesis  and 
some  relevant  theorems  are  established. 

The  thesis  is  concerned  with  the  asymptotic  behavior  of 
sequences  whose  generating  function  has  a  finite  number  of  poles 
or  branch  points  on  the  circle  of  convergence.  A  method  of 
approximating  the  terns  of  such  a  sequence  was  developed  by 
M.  G.  Darboux  in  1878.  His  proof  depended  on  the  order  of 
magnitude  of  the  terms  in  a  trigonometric  series.  Our  proof  of 
the  Darboux  result  is  presented  in  chapter  two.  This  proof  is 
obtained  by  means  of  contour  integration.  It  yields  the  final 
result  in  a  more  convenient  form  and  is,  as  far  as  we  know7, 
unpublished. 

In  spite  of  the  convenience  of  this  method  and  the  large 
number  of  sequences  to  which  it  can  be  applied,  it  does  not  seem 
to  be  well  known.  Some  of  the  examples  given  in  chapter  three 
deal  with  sequences  whose  asymptotic  behavior  has  been  established 
by  some  other  method,  and  to  which  the  method  of  Darboux  is 
applicable.  In  an  of  these  cases,  the  method  used  is  specialized 
and  yields  only  the  dominant  term  of  the  expansion.  On  the  other 
hand,  the  method  of  Darboux  gives  the  complete  asymptotic  behavior 
of  each  sequence  easily. 
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CHAPTER  I 

DEFINITIONS  OF  ASYMPTOTIC  EXPANSIONS 

The  idea  of  an  asymptotic  expansion  of  a  function  began 
with  Poincare'"  in  1886  [  1*1  ]  .  His  definition  is  as  follows* 

DEFINITION  1.1  : 


A  series 


A1  A2 

(1.1)  A0  ♦  -  ♦  - 

z  z2 


is  called  an  asymptotic  expansion  of  the  function  f(z)  ,  in  a 
given  range  of  arg  z ,  if,  for  every  fixed  value  of  n  , 


(1.2)  lim  z3 

I  Z  I  — ^  00 


f(z)  -  A  - 


s  0 


The  notation  used  to  denote  such  a  relationship  between  the 
series  and  f(z)  is 


1  *11 


(1*3)  f(z)  ~  A  . 


The  following  properties  of  an  asymptotic  expansion  are  well 
known* 

If  a  function  f(z)  possesses  an  asymptotic  expansion  in 
Poincare‘S *s  sense,  this  expansion  is  unique.  These  unique 
coefficients  can  be  determined  successively  by  the  equations 


(l.U)  A0  -  lim  f(z) 
l  z  |  —  «> 


11m  z  I  f(z)  -  A 


|Z| 


(1.5) 


. 


.  - 


-  2  - 


(1.6)  ^  =  lim  zn  f(z)  -  Aq  - 

I S I  -*  00  '  z 


Ai 


If 


(1.7)  f(z)-  2  A  z 

r-o 


and 


(1.8)  g(z)-v  2  B  z~r 
r=o 


for  the  same  range  of  values  of  arg  z  ,  then 


(1.9) 


a  f(z)  +  p  g(z)  ~  2  (a  A  +  p  B  )  z”r 
r=o  r  r 


and 


(1.10)  f(z)  g(z)  ~  2  Cr  z"r  , 

r-o 


'where 


(1.11)  C  s  A  B  +  A  B  _+...  +  A  B  . 
r  or  1  r-1  r  o 


If 


(1.12)  f(z)  2  A  z~r  ,  a  <  arg  z  <  p 
r-o  r 


then,  f(z)  -  A  -  — 
o  z 


is  integrable,  and 


(1.13)  F(z) 


Ai  ) 

f(t)  -  Ao  -  “  (  dt 


A_  A  A, 

2  3  It 

~  —  +  — P  +  — ,  +  . . .  ,  a  < 

z  2z  3zJ 


-  A^-l  1 

zn-l  5 


arg  z  <  p 


-  3  - 


Let  f(z)  be  a  differentiable  function.  Then  if 

00 

(l.lU)  f(z)  ~  2  Ar  z“r  ,  a  <  arg  z  <  p 

r=o 

and  if  f*(z)  possesses  an  asymptotic  expansion,  then 

A  2A 

(1.15)  f’(z)  ^  -  — r—  «...  ,  a  <  arg  z  <  p  . 

Z  7? 

It  should  be  noted  that  nothing  has  been  said  about  the 

CO 

convergence  of  the  asymptotic  series  2  A  z~r  ,  From 

r-o  r 

the  point  of  view  of  usefulness,  the  convergence  or  divergence  is 
unimportant . 

If  a  function  f(z)  does  possess  an  asymptotic  expansion, 
this  expansion  may  be  used  to  compute  the  values  of  f(z)  for 
large  values  of  z  ,  the  error  after  n  terms  being  o(l/zn~“). 
In  fact,  in  most  cases,  the  first  few  terms  are  sufficient  for 
computational  purposes,  the  accuracy  increasing  as  |z  |  increases. 

It  is  easily  shown  that  there  exist  functions  that  do  not 
possess  asymptotic  expansions  in  the  Poincare^  sense.  For  example 

it  --if 

(1.16)  zZ/(z  +  1)  =  1  -  z~*  +  z  -  z~Z  +  ... 

(1.17)  (z  +  log  z)"1  =  z-1  |  1  -  (log  z)/z  +  (log  z)2/z2 

expansion  dll 

In  the  first  the  powers  of  z  are  not  integers;  in  the  second, 

A  A 

the  factor  z“^  multiplies  the  entire  series  and  the  numerators 
of  the  terms  of  the  series  are  functions  of  z  rather  than 
constants.  Yet  each  of  these  series  has  many  of  the  desirable 
properties  of  an  asymptotic  expansion. 


? 


■f 


-  h  - 


The  usefulness  of  asymptotic  expansions  has  resulted  in 
numerous  generalizations  of  Poincare’s  concept# 

’One  of  the  most  recent  and  most  elegant  definitions  of  an 
asymptotic  expansion  is  due  to  Professor  A.  Erdelyi.  This 
definition  has  not  yet  been  published  and  was  disclosed  in  a 
private  conversation.  It  is  dependent  on  the  concept  of  an 
asymptotic  sequence  which  was  introduced  by  Erdelyi  [ 1*2  ]  . 

In  any  discussion  to  follow,  the  complex  variable  z  will 
lie  in  a  region  R  and  zQ  will  be  a  limit  point  of  R,  which 
may  or  may  not  lie  in  R  * 


DEFINITION  1.2  : 

The  finite  or  infinite  sequence  of  functions  j  j^(z)  j 

is  called  an  asymptotic  sequence  as  z  -►  z  in  a  region  R  if 

for  each  n  ,  ^n(z)  is  defined  in  R  and  ^p4^(z)  =  °($n(z)) 

as  z  -»  z  in  R. 
o 

Among  the  examples  which  Erdelyi  gives  of  asymptotic  sequences 


are  the  following: 

(i)  (  Z"  Xn  i 

,  where  R1  A  -  >  R1  A. 

3  n+i  n 

for  each  n  , 

and  z  -*  o®  in  the 

sector  0  <  lz/  <  »  ,  1  arg  z 

1  <  rr  , 

(“nz  «•  A  vj  } 

e  z  j  ,  under  the  same  conditions  as  above, 

I  I  IT 

/  arg  z  /  <  -  . 

Two  sequences,  j  ^(z )  j  ,  j  ^(z )  j  such  that 
/n(z)  =  0(  /n(z))  */n(z)  =  °(/n(z))  $  are  said  to  be 

equivalent.  It  is  easily  shown  that  if  j"  j^(z)  }  and  ^  z)  / 


(ii) 

except 


„ 

„ 


, 


n 


. 


1 


are  equivalent  sequences  and  j  /n(z)  f  is  an  asymptotic  sequence, 
then  i  y^n(z)  J  is  also  an  asymptotic  sequence* 

Erdelyi  also  discu  es  parious  methods  of  obtaining,  from 
given  asymptotic  sequences,  new  such  sequences.  For  example, 

(i)  Any  subsequence  of  an  asymptotic  sequence  is  an 
asymptotic  sequence, 

(ii)  If  j  ^n(z)  J  is  an  asymptotic  sequence  and  a  >  0  , 
then  |  |^n(z)|  j  is  an  asymptotic  sequence. 

(iii)  If  j  ^n(z)  j  and  j  y^(z)  J  are  asymptotic 
sequences  containing  the  same  number  of  functions,  then 

j  rfn(z)  y^(z)  j  is  an  asymptotic  sequence. 

We  are  now  in  a  position  to  state  Erdelyi’ s  definition. 


DEFINITION  1.3  s 


Le- 


j  ^r(z)  j  be  an  asymptotic  sequence  as  z— >zQ  in  E, 
Then  the  series  Z  y^(z)  is  said  to  be  an  asymptotic 
expansion  to  N  terms  of  f(z)  as  z  zQ  in  R,  if 


(1.18)  f(z)  *  2  ijy  {%)  +  o(/g(z))  as  z -*z0  ,  s  =  1, 


r=l 


An  immediate  consequence  of  this  definition  is 


(1.19)  ^r+1(z) 


»(/{*)) 


as  z  z 


0  ,  r  =  1,  N-l 


We  use  the  notation 


. 


'  . 


■ 


-  6  - 


N 

(1.20)  f(z)  ~  2  (z) 

r-1 

to  denote  this  relationship. 

If  we  let  y^(z)  =  Ar  ^r(z)  ,  where  the  Ar*s  are 

constants,  we  obtain  Erdelyi’s  published  definition. 

If  a  function  f(z)  possesses  an  asymptotic  expansion  in 
terms  of  a  particular  asymptotic  sequence,  this  expansion  is 
unique.  On  the  other  hand,  one  and  the  same  function  may  have 
asymptotic  expansions  involving  two  different  asymptotic  sequences, 
and  the  two  sequences  need  not  be  equivalent.  An  asymptotic 
expansion  does  not  determine  the  corresponding  function  uniquely, 
since  several  functions  may  possess  the  same  asymptotic  expansion. 
Therefore,  it  seems  desirable  to  place  such  functions  into  the 
same  equivalence  class.  Erdelvi  accomplishes  this  by  the  following 
definition  of  asymptotic  equality  [  1.2  ] 

DEFINITION  l.U  : 

Let  f(z)  and  g(z)  be  two  functions  defined  in  R  and  let 
/rf  (a)  J  be  an  asymptotic  sequence  for  z  -*  z0  in  R.  Then 
f(z)  and  g(z)  are  said  to  be  asymptotically  equal  with  respect 
to  {  f!r(z )/  if 

(1.21)  f(z)  -  g(z)  =  o(^r(z))  as  z  ^z0  in  R,  for  each  r. 

Clearly,  if  f(z)  ~  2  Ar  /r(z)  as  z  ~^z0  and  if  f  (z) 

is  asymptotically  equal  to  g[z)  with  respect  to  \  $zfr(z)  /  , 

then  g (z )  2  Ar  tfr(z)  as  z  zQ  . 


. 


p 


% 
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The  relationship  of  asymptotic  equality  between  functions  is 
obviously  a  true  equivalence  relation. 

Thus  each  asymptotic  sequence  separates  all  functions  into 
classes  of  asymptotically  equal  functions  in  such  a  way  that  all 
members  of  a  particular  class  possess  identical  asymptotic 
expansions. 

The  dependence  of  these  asymptotic  expansions  on  asymptotic 
sequences  can  be  inconvenient.  For  example,  if  f(z)  ~  2  j  , 

the  identical  series 

(1.22)  A«^  +  0  +  +  0  +  A^  +  0  +  •  •  •  , 

is  no  longer  an  asymptotic  expansion  of  f(z)  since  the  sequence 

(1.23)  ^2  *  0  y  ^2  9  0  y  ^  J  0  ,  ••• 

is  not  an  asymptotic  sequence.  Since  this  objection  is  overcome 
by  Erdelyifs  more  general  definition,  we  will  use  the  following 
specialization  of  this  latter  definition  for  the  remainder  of  this 
thesis. 

Let  |Ar(z)  }  be  a  sequence  of  functions  defined  in  R  ,  and 
let  U(z)  be  a  function  such  that  ) U (z ) |  as  z-^z0  in  R. 

Clearly,  f  (U(z))~r  j  is  an  asymptotic  sequence.  The  formal 

CO 

sum  2  Ar(z)/(U(z))r  is  said  to  be  an  asymptotic  expansion 

r-o 

of  f (z)  as  z  -*  z0  if,  for  all  n  , 

(1.2U)  f(8)  -  Z  Ar(z)/(u(z))r  =  0(l/[u(z)]a)  as  z  -  z0  . 

r=o 


* 


e. 


-  3  - 


Again  we  use  the  notation 


(1.25)-  f(z)~  2  Ar(Z)/(U(Z))- 

r=*o 


The  usual  meaning  is  attached  to 


(1.26)  f(z)~g(z)  2  Ar(z)/(n(z))' 

r-o 


THEOREM  1.1 


If  f(z)  ~  2  A  (z)/(U(z))r  as  z  -5-z0  in  R, 


r-o 


then  lim  An(z)/u(z)  =  0  ,  for  all  n  „ 

z 


PROOF: 


As  a  result  of  our  definition,  we  have 


(1.27) 

f(z) 

-  2 

r*o 

and 

(1.28) 

f(z) 

n+1 

-  2 

r=o 

Hence, 

(1.29) 

A  .  (z 
n+1' 

)/(U(z 

and 

0(1/  [tf(z)]n) 


\]  n+1- 


n+1 
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(1.30)  An+1(z)/U(z)  =  o(l)  , 

for  all  n  ,  i.e., 

(1.31)  lim  An(z)/U(z)  =  0 

2  -*  Z0 

For  convenience  we  shall  adopt  the  following  specialization 
of  the  Frdelyi  definition  of  asymptotic  equivalence. 

Two  functions  f(z)  and  g(z)  will  be  called  asymptotically 
equivalent  with  respect  to  U(z)  at  z  =  zQ  ,  if  positive 
constants  a  and  b  ,  and  a  neighborhood  N  of  zQ  exist  such 
that,  for  all  z  common  to  N  and  R  , 

(1.32)  I  f(z)  -  g(z)  |  <  a  e_u'  1  . 

1 

-  z 

For  example,  the  functions  f(z)  -  z  and  g(z)  -  z  +  z 
are  asymptotically  equivalent  with  respect  to  U(z)  =  In  z  at 
z  =  .  However,  the  same  two  functions  are  not  asymptotically 

equivalent  with  respect  to  U(z)  =  z  at  z  =  «•  . 

This  again  is  obviously  an  equivalence  relation. 

THEOREM  1.2  : 

A  necessary  and  sufficient  condition  for  two  functions 
f(z)  and  g(z)  to  be  asymptotically  equivalent  with  respect  to 
U(z)  at  zq  ,  is  that  there  exist  a  sequence  of  positive  numbers 
j an  ]  9  aAnumber  ncM,  and  a  fixed  neighborhood  of  zQ  for  which 

positive 

00 

(i)  Z  an  cn  converges 

n=o 

(ii)  I U | n  |  f(z)  -  g(z)  |  <  an  nl  for  all  z  in  the 

intersection  of  R 
and  the  neighborhood. 


' 


PROOF: 


Let  us  first  prove  the  necessity*  There  exist  constants  a 
and  b  ,  and  a  neighborhood  N  of  z0  such  that,  for  z  common 
to  N  and  R, 

(1.33)  |  f(z)  -  g(z)  |  <  a  e"blUl  . 

Therefore , 

(1.3U)  |U|n  |  f(z)  -  g(z)|  <  a  |U|n  e“bl,jl 

If  we  consider,  for  a  moment,  the  function  w  «  x31  e“DX  , 
x  >  0  ,  we  see  that  since 

(1.35)  £  =  xn"1  e'bX  (n  -  bx)  » 

w  must  have  a  horizontal  tangent  at  x  -  n/b  ,  This  point 
obviously  yields  a  maximum  since  w  is  a  positive,  continuous 
function  which  tends  to  zero  as  x  tends  to  °°  or  0  * 

Applying  this,  we  obtain 

(1.36)  |u|n  |  f(z)  -  g(z)  j  < 


If  we  choose  a  =  . -n 

rr  n  b 

theorem  is  satisfied*  By  a  simple  application  of  the  ratio  test. 


ef 


a  e~n 

n*  bn 


\/T  it  n  b1 


ni 


,  condition  (ii)  of  the 


11  - 


(1.37) 


Z  ^  cJ 
n-1 


VSf  2 


cn 
t~Z  n 


n-1  n 


is  seen  to  be  convergent  for  c  <  b  .  Thus  condition  (i)  is  also 
satisfied. 

To  prove  the  sufficiency,  we  deduce  from 


(1.38)  |u|n  |  f(z)  -  g(z)  \  <^111 


that 

(1.39)  2  'S-lh  |  f(a)  .  g(z)  |  .  2  an  cn  a 

n=o  n  ’  n=o 

and 

(l.Uo)  |  f(zj  -  g(z)  |  <  a  e-c 1  1  . 

Thus,  from  the  definition  of  asymptotic  equivalence,  it  follows 
that  f(z)  and  g(z)  are  asymptotically  equivalent  with  respect 
to  U  at  z 

o 

The  following  two  theorems  connect  the  concept  of  asymptotic 
equivalence  to  that  of  asymptotic  expansions. 


THEOREM  1.3  : 

If  the  functions  f(z)  and  g(z)  are  asymptotically 
equivalent  with  respect  to  U(z)  at  z  *  zq  ,  then 

f(z)  g(z)  as  Z  -»  z0  . 


PROOF: 

This  follows  trivially  since  there  exist  constants  a  and 
b  and  a  neighborhood  N  of  zQ  such  that,  for  all  z  common  to 


N  and  R 


> 


4 
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(l*Ul)  |  f(z)  -  g(z)  |  a  e~k|ul  «  o(l/jjn)  for  all  n,  as  z 


THEOREM  l.U  : 

00 

If  f(z)  ~  Z  A  (z)/(lT(z))r  as  z  -*  zQ  ,  and  if  the 
r=o 

corresponding  terms  of  the  two  sequences  |Ar(z)|  and  jB_^(z)J 
are  asymptotically  equivalent  with  respect  to  U(z)  at  z  =  zQ  , 

CO 

then  f(z)  ~  Z  3  (z)/(U(z))r  as  z  -*>  zQ  „ 

r^o 

PROOF: 

As  a  result  of  the  preceding  theorem, 

(1.1*2)  Aq(z)  -  Bn(z)  =  o(l/un)  for  all  n,  as  z  zQ  * 

Consequently, 

f(z)  -  2  Br(z)/(U(Z))r 

r=o 

=  f(z)  -  2  {  Ar(z)  +  [Br(z)  -  Ar(z)]  ]  /(U(z))r 

r=o  ' 

(1.U3) 

=  f(z)  -  2  Ar(z)/(u(z))r  +  0(1/11°) 

r=o 

=  o(l/Un)  as  z  z0  . 

Therefore, 

00 

(1.1*U)  f(z)  A.  2  Br(z)/(u(z))r  as  z zQ  . 

r=o 


- 


' 
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chapter  II 


THE  METHOD  OF  DARBOUX 

In  the  year  1878  Darboux  considered  the  problem  of 
determining  the  asymptotic  behavior  of  a  sequence  of  numbers,  [* 
In  our  present  chapter  we  shall  give  an  alternate  derivation 
of  his  result  under  slightly  different  assumptions. 

We  shall  assume  that  a  generating  function  f(z)  exists 
for  the  sequence  j  an  }  of  the  form 

CD 

(2.1)  f(z)  =  Z  an  zn 

n=o 

and  that  this  generating  function  has  the  following  properties : 


(A)  f(z)  is  a  regular  function  of  z  for  |s|  <  R  , 

(B)  f(z)  has  a  finite  number  of  singularities  on  its  circle 

of  convergence  at  z  ~  a.  ,  j  ~  l,,,.,m 

J 


(C)  In  the  neighborhood  of  the  singularity  z  =  a.  ,  f(z) 


has  the  form 


-P 


(2.2)  f(z)  =  (z  -  a.)  J  g  (z)  +  h.(z) 

J  J 


where  g  (z)  and  h.(z)  are  regular  functions  of  z  in  a 
j  3 

neighborhood  of  z  »  a.  ,  and  R1  p.  >  0  ,  j  -  l,,,.5m  • 

5  J 

Our  assumptions  imply  that  the  singularities  of  f(z) 
which  lie  on  the  circle  of  convergence  are  either  poles  or 
branch  points. 

By  our  assumptions,  g.(z)  has  a  Taylor  expansion  about 

0 

z  =  a.  of  the  form 
J 


♦ 


. 


0 


, .  *  ■  *  ' n  • 


. 


. 


« 


(2.3) 


-  lU  - 

g.(z)  -  2  (o  )  (z  -  a  )r  /rl  , 

2  r=o  3  ^  3 

which  We  write  as 

k-1  v 

( 2 *U)  g.(z)  *  Z  (a  )(z  -  a  )  /rl  +  (z  -  a.)k  K  (z) 

«J  rso  j  J  J  J  J 

Let  s  =  max  (Rl  p.)  ,  j  =  l,.,.,m  .  Since  k  is 
J 

arbitrary,  we  may  choose  k  >  [s]  ,  where  [s]  denotes  the 

largest  integer  contained  in  s,  Further,  we  define  an 
integer  q  by 

(2.5)  q  s  k  -  [s] 

From  (2.U),  we  can  write 

k-1 

(2.6)  f(z)  =  2  (a  )(z  -  a.)r_Pj  /rl 

r=o  J  J  J 

k-p. 

+  (z  -  a.)  J  K  (z)  +  h.(z)  , 

J  5  3 

in  a  neighborhood  of  z  *  a.  * 

0 

If  we  make  appropriate  cuts  in  the  z  plane,  it  is 
possible  to  find  a  region  of  the  shape  shown  below,  within 
which  f(z)  is  regular*  This  region  shall  be  denoted  by 
U(a  ,  S  ),  or  simply  U  *  If  desired,  U  could  be  made  to 
extend  out  to,  but  not  including,  the  next  singularity  of 
f(z).  However,  in  our  proof,  both  ”aM  and  M  £"  will 
eventually  be  made  to  approach  zero. 


-  .. 
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The  expression  D^(z),  given  by 


(2.7)  Dk(z) 


m  k-1  ,  , 

2  2  g\r'  (a.)  (z  -  a.)  “Pj  /r , 

j=l  r=o  J  J  J 


is  called  the  Darboux  approximant  of  f(z).  Clearly,  Dk(z) 
is  regular  within  the  region  U,  and  has  a  Maclaurin  expansion 
convergent  within  |z|  <  R.  Thus 

00 

(2.8)  D^(z)  s  2  zn  . 

n=o 

The  function  F(z),  defined  by 


~r 


. 

* 


‘ 
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(2.9)  F(z)  -  f(z)  -  Dk(z) 


=  Z 
n=o 


(a  -  a  .  ) 
v  n  nky 


9 


is  also  regular  within  U.  In  a  neighborhood  of  z  =  a-  , 

J 

F(z)  has  the  form 


k-1 

(2.10)  F(z)  =  Z  g[T\a  )(z,  »  a  )r~pj  /rl  +  (z  -  ct.)k~Pj  K.(z) 
r=o  J  J  J  3  3 

m  k-1  j  .  ^  _ 

+  h  (z)  -  Z  Z  g .  (a.) (z  -  a  )  Fj  /rl 

J  j~l  r~o  3  ° 

=  (z  -  a  )k'Pj  K  (z)  *  H  (z)  , 

J  3 


where  H.(z)  is  regular  in  this  neighborhood. 
0 


It  thus  follows  that 


(2.11)  F(q)(z)  =  (z  -  a,)k"pj"q  K.  (z)  +  H<q)(z) 

j  39  q  3 

=  (z  -  a.f~pj  K.  (z)  +  H^q)(z) 

3  3 

where  the  functions  K.  (z)  are  regular.  Thus  F(z)  has 

0 

q-1  bounded  derivatives  in  a  neighborhood  of  z  =  a.  . 

3 

From  (2.9),  we  have 


(2.12)  an  -  &nk  =  (2  tt  I)'1  j*  F(z)  z'(ml)  dz  , 

c 

where  C  is  any  contour  within  the  domain  of  regularity  of 
F(z)  .  Integration  by  parts  q  times  yields 
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M)  an-ank  “ 


X (n-q)l  f  F^(z) 

tt  i  nl  J  n-q+1 


dz 


Let  us  divide  the  contour  C  up  into  three  parts? 


I  The  small  circles  enclosing  the  singularities. 


II  The  radial  lines 


III  The  outer  circle  of  radius  R  +  $ 


We  shall  now  choose  a  «  R/n  *  and  £  =  R/  \/n  . 


On  the  small  circle  enclosing  the  singularity  z  =  a.  * 

3 


we  have 


(2.lU)  z  «  a.  +  a  e 
J 


ie 


Therefore*  from  (2.11)  constants  M  and  N  exist  such  that 


(2.15) 


v(q) 


z) 


W  -R1  PJ  M  + 


and 


(2.16) 


,(q) 


(z) 


n-q+1 


dz 


[s]  -R1  p  .+1 
a  3  m  +  a  N 


(R  -  a) 


n-q+1 


2rr 


Rn 


On  the  radial  lines* 


(2.18)  z  =  reiargaj 


Again  from  (2.11)*  we  have 


(2.19) 


v(q) 


00 


(r  -  R) 


[s]-Rl  p. 


M  +  N 

i  i 
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Thus, 


(2.20)- 


fW(z) 


R  *  £ 


n-q-1 


dz 


II 


(r  -  R) 


[s]-Rl 


\  +  N1 


,n-q+l 


dr 


R+a 


a 


[s]-Rl  p.+l  W-Rl  p.+l  f 


)  \  (S  -  a) 


n-q+1 


On  the  large  circle  of  radius  R  +  £  , 

(2.21)  z  =  (R  +  S  )  eie  =  R(1  +  Jj»  )  eie 


Thus, 


(2.22)  z  -  a.  =  R  f  (1  +  F  )  eie  -  e1  arg 


1 
Vn 


and 


(2,23)  vff  -  I  z  ‘  aj  I  -  K(2  +  y5  } 


If  [s]  -  R1  p.  >  0  , 

0 


(2.2k)  |  ,<«>(,)  |  <  I E(2  *  i)|I’1‘F1P3 


V® 


M2  +  N2  , 


and 


(2.25) 


I 


III 


F^(z) 

n-q+1 


dz 


R(2  +  =  ) 

<  v. 


^  j  [s]-Rl 


*  Mp  +  N 
2  2 


Rn_<3  (1  +  i  )“"q 

vn 


2tr 


o 
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If  [s]  -  Rl  p.  <  0  5 

J 


(2.26)  ' 


and 


F^(z) 


-  \Jn 


[s]  -  RI  p 


M2  *  N2  , 


(2.27) 


F(c^(z) 


dz 


,n-q+l 


III 


R\[s  ]  -  Rl  p.  +  1 


2  Mo  +  F  N. 


y/n 


^n-q+1  (  1  (  x  +  .  1  )n-q  J 


2tt 


Vn 


These  results  and  (2.13)  show  that 
(2.28)  • 


(n-q )  l  q  / 1 


Using  (2.7)  and  (2.8)  ,  we  obtain 


(2.29)  a. 


nk 


i 

nj 


dz 


Bk(2) 


m  k-1 


•Wfcx  1  f-a  )r-P.rn  r  (n+p^-r) 


nj 


2  2  g)  (O  (-oj  *3 

3=1  r=o 


r.'  r(p  -r) 


Thus 


(2.30)  s  -  -- 


(-Dn  ;  k-’- 


j=l  r=o 


i  gSr)(a.)(-a.)r-P3-n  r(n+p j-r)  ,  ? 

rl  P*(p^-r)  nl 


This  provides  us  with  an  asymptotic  expansion  for  an  3  since 
for  each  k  >  [s]  3 


(2. .31)  an  =  (_Dn  2  k2X  gjr)(a3)(~C3)r~Pj  / R  P^) 


P(n+s) 


+  <s=a2l  .  o  (i) 

T(n+s) 


3=1  r=0 


r  J  r(p,-r)  V  3/  r(n+s) 


! 


1 


I  ■ 


1 


-  20  - 


which  can  be  put  into  the  form 


(2.32) 


nt  Rn 

r(n+s) 


A0(n)  + 


Ai(n) 

n 


Ak-l(n) 


n- 


k-1 


Hence 


(2.33) 


(-l)n 


m 
2 

nJ  i=i 


2 

r=o 


(a.x-cor-prn 


r(n+p.-r) 


rJ  r(Pj' 


•r) 


In  the  event  that  is  a  pole  of  order  p  ,  it  is 

easily  seen  that  its  contribution  to  a^  and,  therefore  also 
to  the  asymptotic  expansion  of  aR  ,  contains  only  p^ 
terms , 


-  21  - 

CHAPTER  III 


APPLICATIONS 


As  a  first  example,  we  shall  obtain  the  asymptotic 
expansion  of  Legendre  polynomials  of  large  degree*  The 
generating  function  of  the  Legendre  polynomials  is 

i  °0 

(3.1)  I  1  -  2tz  +  z2  I  "  2  =  2  P  (t)  zn  . 

n=o 

We  shall  first  choose  t  =  cosh  x  ,  where  x  is  real  and 
positive*  Thus,  (3*1)  becomes 

(3.2)  2  Pn(cosh  x)  zn  =  j  1  »  (ex  +  e“x)  z  +  z2  f 

n=c  * 


The  minus  sign  comes  from  choosing  the  branches  of  the 

i  i 

functions  j  z  -  e“^  j  ~  2  and  |  z  -  ex  ]  ”  z  that  are 
Ex  «» 

equal  to  -i  e  and  -i  e  '  respectively  at  z  =  0  * 

Clearly,  the  first  singularity  is  at  z  =  e”x  *  Therefore, 
in  the  notation  of  the  preceding  chapter. 


i 


(3.3) 

g(z)  =  -f 

*  -  -  1  , 

(3.U) 

g(r)(Z)  « 

r(r+l)  ,  .r+l  , 

■p/i^  (*“1)  (z  —  e  )  , 

and 

(3.5) 

g(r) (e~X)  - 

r  (r+l)  ,  .  >r+l  .  „  .  .  ,-fr+i) 

■p  T,.  (-1)  (  -  2  sinh  x  ) 

1  i  2) 

From  (2.33),  we  have 


N\'-k 


f 
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(i+n)  :  r(rx-l)  r  (n-r+J) 


(3.6)  P  (cosh  x),,  a__  Z  A  A'-*;.  .f.  -V2=-V*;  (_].)*  e~™  (2  sinh  xp"*. 
nlVd)  r=o  T(l-r)  rl 


,-(m 


As  a  result  of  the  identity  V  (w)  F  (1-w)  *  Tr/sin(iT  w)  , 


(3.7)  r  (i  -r)  r  (I  +r)  =  rr(-l)r  ,  and  T  (n-r+|)  T  (r-n+J)  =  tr(-l) 


n-r 


Thus 

(3.8)  Pn(cosh  x).  S^L  (=^  2 

V2rr  sinh  x  nj  r=G  rj  r  (r-n+|) 


-e 


-x 


2  sinh  x 


x(|+n) 


(-Dn  ril)  T(i)  „ / j  i 


V2tt  sinh  x  n.f  P(-n+|) 
r(n+|) 

nl  v^tt  sinh  x 


F(i»b-“*25  -e“X/(2  sinh  x)) 


F(  1  >  1 5  -n+  {  ; 


-X 


2  sinh  x 


If  we  now  choose  t  « 

cos  6 

,  (3*1)  becomes 

eo 

(3.9)  2  Pn(cos  e)  zn 

=  1 

1  -  z(e^@  +  e"iQ)  +  z2  } 

i 

n=o 

=  -{ 

z  -  J  *  |z  -  e"j  e  \ 

i 

% 

Again,  the  minus  sign  comes 

from  choosing  the  branches 

of 

the  functions  \ z  -  e1^  ( 

i 

-  z 

and  /  z  -  ^  "  2 

that 

-  -  i9 

are  equal  to  -i  e  2 

and 

|i© 

-i  e  respectively  at  z 

i  ft 

In  this  case,  there  are  two  singularities,  ax  s  e  and 

*""1  Q 

a2  ~  e  ,  both  on  the  unit  circle.  Since  the  contribution 
of  one  singularity  will  clearly  be  the  complex  conjugate  of  the 
contribution  of  the  other,  we  need  only  consider 


t. 
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i9 

ai  e  e 

5  gi(z)  “ 

-  {  z  -  e-'9  \  ~  *  . 

Then 

(3.10)  ' 

gir)(z)  = 

,  ,vr+l  r  (r+i)  c 

Hi)  ‘ 

z  .  e-xe  ^  -(r+i)  f 

and 

(3.11) 

g^r)  (ei9) 

/  n %r+l  r (r+i) 

r(t) 

[  21  sin  6  | 

Therefore,  from  (2.33)?  we  have 


(-l)n  09  r(n+£-r)r(r+i)  ,  Xr+1,  3©H‘n,  ;(r+i) 

(3.12)  Pn(oos  e)~  2  RX  ^  r,  r(f)"r~^-r')(  )  (_6  (asin9) 


.  tf  J  tfcMiiijSlL HHeh, 

n*  r=o  r (r-n+ 1)  rjrm 


•[©(r-i-n)+r(l-2r)] 


(ii  )nfZH 

nl  /tt  sin  © 


^jrCr+i)}2  cos  j  9(n-r+  -f)  +  jj  (2r-l )j 

r”°  F(r-n+|)  rl  (2  sin  ©)r 


) 

■  :  i 


i 
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As  a  second  example  we  shall  discuss  the  asymptotic  behavior 


of  the  sequence  j  uR  j  whose  generating  function  is  given  by 


x  + 


3.13) 


(3.1U) 


i  x2 

00 

x 

' 

2 

. 

n=o  nl 

% 

=  b 

n 

/  2n  ,  we 

have 

x^ 

ex  + 

2  x2 

2n 

nl 

V  i  - 

X 

i  e 


x  +  £X‘ 


\/~xT-  1 


The  generating  function  has  a  branch  point  at  x  -  1  . 
the  notation  of  chapter  two. 


In 


(3.15)  g(x)  »  i  e 


X  +  2  X‘ 


and 


(3.16)  g(r)(l) 


~  (i  © 

dx2 


x  +  \  x2  > 


J  x=l 


dr  .  (x+1)  +  2  (x+1) 2 

—  i  e 


dxr 


A— w 


-  i  e 


dxr 


d*  i  x2  +  2x 


x=0 


This  derivative  can  be  evaluated  by  the  use  of  the  Hermite 
polynomials.  The  generating  function  of  these  polynomials  is 


2tz  -  z2 


2  H  (t) 

n  nj 


n-o 


(3.17)  e 


' 


, 


If  z  *  i  x  /  ^  ,  t  ■  -  i 


9 


(3.17)  becomes 


(3.13)  • 


|x2  +  2x 


2  Hn(-  y/2  i) 
n=o 


xn 

n  l 


Thus 


(3.19) 


e  2  x2  +  2x 

dx* 


x=0 


and 

(3.20)  g(r)(D 


Z  /  i\  r 

le  (v?  Hr(-V2i) 


The  Hermite  polynomials  are  given  explicitly  by 


(3.21)  Hjx) 


[I] 


nl  2 
m=o 


(-l)m  (2x)n-2w 

ml  (n-2m) 1 


Hence 


(3.22)  Hr(-  v^  i) 


[f]  44^  . 

ft—  o  m  •  ( r— 2m )  1 


and 


(3.23)  H  (-  y/S-  i) 


y-  yfl  i) 


1  , 

•  2  V?  i  , 


H2(-  i)  -  -  10  , 


H,(-  y/S  i) 


28  yjz  i  . 


There  fore ,  from  (2.33) 


(3.2li)  fc  ~  2n  e1  Z  L  i]  —  ”~r+  ~  H  (-  y/2  i) 
n  r=0  \  yfz)  p(t-r)  ri  ^  ^  ^ 


If  we  compute  b^  and  b 2 ^  using  the  first  four  terms  of 


this  expansion,  i.e., 


I 


f 
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(3.25) 

b  ^  2n  e2 
n 

r(m  1) 
T(  £ ) 

l1 

we  obtain 

(3.26) 

bl5 

2.997 

x  10: 

b25 

2.7352 

x  id 

1? 


2n-l 


2(2n-l)(2n-3) 

35 


(2n-l)(2n-3)(2n-5) 


+  . . . 


16 


32 


To  determine  the  accuracy  of  this  expansion,  let  us  now 
establish  a  recurrence  relationship  for  the  sequence  j  bn  j 
and  find  the  exact  values  of  these  coefficients. 

Differentiation  of  (3 *13)  with  respect  to  x  yields 


(3.27) 


x  +  4  x2 

e  * 


v  1  -  x 


(1  +  x)  + 


1 

X  +  gx2 

e 

2(1  -  x) 


CO 

z 

n-o 


uml 
n  J 


x31 


which  simplifies  to 


(3.28) 


Z 

n~c 


nl 


Z 

n-o 


x11  (i  -  X)  . 

ni 


Hence  we  have  the  recurrence  relationship 

(3.29)  (  I  +  n)  un  -  n  (n-1)  un_2  =  un+1  , 
and  the  related  equation 

(3.30)  bn+1  =  (?n  +  3)  bn  -  8  n(n-l)  bn_2 

By  expanding  the  left  hand  side  of  (3.13),  we  obtain  the 


initial  values 


5 


I. 


c. 


f 


-  27  - 

3  15 

(3.31)  u0  =  l,u1  =  2>u2  =  'iT  » 
and 


(3.32)  bQ  -  1  ,  \  -  3  ,  b2  »  15 


Successive  use  of 

(3*30),  together  with  (3.32),  yield; 

the 

values  given  below. 

b^ 

1 

0 

bi 

TS 

3 

b2 

- 

15 

b3 

= 

89 

bU 

= 

657 

b5 

“ 

5  787 

b6 

53 

60  991 

fc7 

as 

757  185 

b8 

= 

10  927  713 

b9 

= 

180  302  579 

b10 

SB 

N 

3  350  213  399 

bll 

S 

69  18?  005  U17 

b12 

= 

1  571  008  865  905 

b13 

- 

38  879  lai  706  891 

biU 

= 

i 

Ola,  157  556  739  U23 

bl5 

= 

29 

988  U95  350  16U  U33 

b16 

=r 

92U 

302  93h  887  8U9  U09 

b17 

3 

30  351 

580  212  135  037  155 

b18 

S 

l  057  753 

501  967  038  568  527 

b19 

s 

38  989  692 

992  109  0H8  819  321 
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b20  = 

l  515 

b21 

61  952 

b22  = 

2  656  855 

b23  = 

119  270  771 

b2U  = 

5  593  U8U  272 

b25  = 

273  535  025  652 

535  U89  2X6  06 9  539  936  081 
U55  390  311  192  968  929  U03 
I2li  110  517  279  568  90U  575 
665  051  719  120  13U  759  6h9 
167  55U  527  7U8  376  999  U57 
U73  236  608  590  9h 4  369  107 


A  comparison  of  the  results  shows  that  is  determined 

to  within  a  0*06  percentage  error,  and  b2^  to  within  a  0*006 
percentage  error. 


m 
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Example  three  is  concerned  with  the  following  problem:  If 

there  be  n  straight  lines  in  one  plane,  no  three  of  which  meet 
in  a  point,  what  is  the  number  of  groups  of  n  of  their  points 
of  intersection,  in  each  of  which  no  three  points  lie  in  one  of 
the  straight  lines? 

Let  us  denote  the  number  of  such  groups  by  .  Robin 

established  the  following  recurrence  relationship 

=  ng  +  (o)  g  , 

n  12/  n-2 

,  g  =  0  ,  g  =  1 

1  2  3 

He  introduced  u  ,  where  g  -  |  (n-l)l  u  ,  and 

n  n  n 

obtained  the  corresponding  relationship 


Robinson 

(3.35) 

g 

n+1 

where 

g_  : 

(3.36)  vi. 


u. 


n+3 


n+2 


n 

2n 


with  u  =  u  =  0  , 

1  2  ’ 

u  = 

3 

\ 

s 

U5  =  1 

• 

Robinson  proved  that  the 

number 

b 

= 

iim  u2 
n  n 

/n 

existed  and  consequently,  that 

^n  - 

/ 

tt  b 

2 

nn  e“n 

.  By 

using  the  sequences  \  u2 

1  n+1 

-  u2  } 
n  ) 

and 

1  Un  /(n  ' 

■M 

each  of  which  converges  to  b,  he  was  able  to  establish 
b  *  0.28U,098, . . .  ,  and  expressed  an  interest  in  finding  an 

exact  expression  for  b  . 

The  May,  19^2  issue  of  the  American  Mathematical  Monthly 
discussed  the  consequent  articles  that  had  been  submitted. 

Various  methods  had  been  used  to  establish  that  Robinson's 

! 

constant  b  is  given  by  b  =  U  e  /  rr  .  The  most  common 
treatment  involved  the  generating  function 


« 


. 

' 


-  30  - 


(3.37)  f(x)  -  2  2  —  x11 

n=3  nl 


Z 

n=3 


n 


x11 


An  explicit  expression  for  f(x)  was  found  from  the 
recurrence  relationship  (3.36)  to  be 

1  O  i 

-  4  .  -  2  x 

(3.38)  f(x)  =  "  ^-:^=  -  2 

V  1  -  X 


Three  of  the  authors  then  used  the  following  theorem:  If 

OO  CO 

f(x)  =  Z  an  x11  ,  |x  I  <  1  ,  and  g(x)  =  Z  bn  x11  , 

n=o  n=o 

OO 

I  x  I  <  1  ,  where  Z  bn  diverges,  and  if  lim  a^  /  bn 

n-o  n-*  00 

then  lim  f(x)/g(x)  =  C  * 

x  — >  1 

Applying  this  to  the  generating  function  f(x),  with 

“  2,  f(x)  __ 

g(x)  s  (1-x)  ,  they  established  that  lim  — — -  i  *  \jv  "b 

x  1  C1"*)  5 

3 

But  from  (3*38)  this  limit  also  equals  2  e *  ,  whence 

.?  3 

b  =  U  e  2  /tt  ,  and  g  ^  ^2  nn  e~n~  * 

One  treatment  employed  contour  integration  and  yielded,  in 
addition  to  the  exact  value  for  b  ,  an  explicit  expression  for  un  • 
Let  us  establish  a  more  complete  asymptotic  expansion  for 
g^  using  the  method  of  Barboux.  We  will  also  make  use  of  the 
generating  function 

(3.39)  f(x)  =  2  z  —  x* 

n=3 

By  using  (3.35),  it  is  easily  shown  that 


(3.1:0)  (l-x)f'(x) 


l  *2  f(x) 


1 


f  -  • 


■ 


i 


. 


I 


-  31 


Since  f(0)  =  0  ,  (3.1+0)  yields 

-  |x2  -  |x 

(3J4I)  '  f(x)  =  2  i 


\Tx  -1 


-  2 


The  function  f (x)  has  a  branch  point  at  x  =  1  .  In  the 
notation  of  chapter  two. 


(3.1*2)  g(x)  =  2  i  e 


••  £x2  -  g x 


and 


(3.U3)  g(r)(l)  - 


2  i 


2  i 


dr  ~  -  gx 

—  e 


L  dx 


dx1 


2  i  el  * 


x=l 

-  |(x+l)2  -  f  (x+l) 


x=0 


U2 


Xfc  -  X 


dxr 


x=Q 


If  z  *  |  x  3  t  -  -  1  ,  (3.17)  becomes 


(3.U h) 


-  Z  X2  -  X 


n 


n-o  2iA  nj 


Thus 


(3.W) 


dx 


(e-^2-b 


Hr(-1) 


-  x=0 


and 


(3.U6) 


g(r)(D 


Hr(-1) 


2  i  e 


Using  (3.21),  we  obtain 


~ 


- 
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(3.U7)  H  (-1) 
r 


rl 


[fi 


^  ^)r”m  2r“^m 

ml  (r-2m) J 


3 


and 


(3.U8)  H  (-1)  =  X 

o 

H^-l)  «  -  2 

H2(-l)  =  2 

Hj(-l)  =  h 


Hence  from  (2.33) 


(3.te,  j  (-  if  $=fi>  v- 


1) 


(3.50) 


8n~  6 


;i  r  (n+  i ) 

rd) 


2n-l  U(2n-l)(2n-3) 


U(2n-l)(2n-3)(2n-5) 


The  dominant  term  of  this  expansion,  together  with  Stirling’s 
formula,  yields  the  result  previously  discussed,  i.e., 

.3 

(3. 5D  gn  ~  e  *  yJz  nn  e~n 

Using  first  (3.51)  ,  then  (3*50)  to  compute  g-M  and  g ^  , 
we  obtain 


(3.52) 

gl5  = 

8.9  x  10 

g25 

8.2  x  10 

10 

23 


and 


-  33  - 


(3. S3)  =  8,625  x  107 

g25  =  806,127  x  1018 

respectively. 

The  exact  values  of  these  coefficients  are  determined  through 
successive  use  of  (3*35)  and  are  given  below. 


gl 

s 

0 

g2 

ss 

0 

g3 

ss 

1 

s 

3 

% 

a 

12 

g6 

= 

70 

g7 

a 

U65 

g8 

s 

3 

507 

g9 

= 

30 

016 

g10 

= 

286 

88)4 

gll 

- 

3 

026 

655 

g12 

as 

3U 

9UU 

085 

g13 

U38 

263 

36U 

glU 

- 

5  933 

502 

822 

gl5 

as 

86  2U8 

951 

2U3 

g16 

- 

1  339  751 

921 

865 

g17 

s 

22  lU8  051 

088 

I48O 

g18 

a 

388  2h6  725 

873 

208 

g19 

as 

7  193  U23  109 

763 

089 

g20 

s 

1U0  U62  355  821 

628 

771 

g21 

2  883  013  99k  3U8 

U8U 

9U0 

-  3U  - 


S22  ■ 

62 

e23 

1  397 

g2U  =  ' 

32  87U 

g25 

806  125 

053  912  73U  368  1*32  U30 
632  88U  350  901  75 9  561 
958  880  6U0  907  159  723 
893  050  067  U59  18U  032 


Comparing  these  results,  we  find  and  g  determined 

respectively  to  within  4*  and  2-  percentage  errors  by  (3*51) 
and  to  within  0.001  and  0.0001  percentage  errors  by  (3*50)  . 

We  shall  now  establish  an  explicit  expression  for  gn  » 
From  (3*3 9)  arid  (3.U1)  , 


(3.5W  g. 


dx* 


1  2  i 
4X2  -  2X 


/"x  1 


-  1 


x=Q 


n 

i  2 
r=o 


(?) 


,n-r 


ip  1 

-  4.  YS  »2X 


dx 


n-r 


x=0 


dx3 


(x  -  1) 


1 

~  2 


x=0 , 


If  z  =  2  x  ,  t  =  -  ^  ,  (3*17)  becomes 


(3.55) 


I  2  1 

-4X2  -  2X 


2  Mil) 

n=0  211  nj 


x11 


Thus 


(3.56) 


■n-r 


1  0 
4X2 


1 


dx 


n-r 


Hn-r(-S) 


>n-r 


x=0 


where,  from  (3.21) 


i. 
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(3.57) 


C-z) 

n-r 


t"]  (_i)n-r-m 


(n-r)  ’  2 

m=o  ml  (n-r-2m)j 


Since 


(3.58) 


dxr 


(x  -  1)_* 


x=0 


i  (2r)» 
22r  rj 


we  have 


(3.59) 


3  S  (V) 

^  r=o 


n-r 


(n-r) J  2J 


n!  n  /2  r 

7  2 

2  r~o 


2  )  i 

\  r  / 


x  r-fi 


(-D 


2r  m=o  mJ  (n_r 


-r-m 

-2m)  l 


I 


Example  four  is  concerned  with  the  number  of  distinct  terms 
in  a  symmetrical  or  partially  symmetrical  determinant. 

In  dealing  with  this  problem.  Professor  A.  Cayley  [s.z] 
represented  the  terms  of  a  determinant  as  duads,  and  the  determinant 
as  a  bicolumn.  For  instance, 


ab  ap f  aq 1 
bb  bp  ’  bq  ’ 
pb  pp 1  pq 1 
qb  qp'  qq'  , 

Thus  a  determinant  whose  bicolumn  contains  such  symbols  as 

pp*  and  qq’  ,rq.(each  letter  p,  q,  ...  being  distinct  from  each 

letter  p*,  q*,  ...  )  is  partially  symmetrical,  while  a 

(  aa  ) 

determinant  such  as  <  bb  |  is  wholly  symmetrical* 

I  cc  J 

Cayley  defined  a  determinant  (m,n)  to  be  one  whose  bicolumn 
had  m  rows  aa,  bb,  ...  and  n  rows  ppT,  qq1,  ...  ,  and 
denoted  the  number  of  distinct  terms  in  the  developed  expression 
of  such  a  determinant  by  <f>  (m,n)  . 

He  established  the  recurrence  relationships 

(3.61)  (m,n)  =  m  0  (m-l,n)  +  n  <f>  (m,n-l)  ?  n  ^  1 

(3.62)  f  (m,0)  =  f  (m-1,0)  +  (m-l)  j>  (m-2,0)  +  (m-l)  (m-2 )  (m-3 ,1) 

and  used  these  equations  together  with  the  initial  values 

i>  (0,0)  =  1,  0(1,0)  =  0(0,1)  =  1  ,  0(2,0)=  0(1,1)  =0(0,2)  =  2, 

to  calculate  values  of  <f>  (m,n)  in  the  following  order 


(3.60) 


represented  the  determinant 


aa 

ba 

pa 

qa 


where  the  duads  had  the  property  that  sr  =  rs  . 


* 

, 


; 


,  ' 


t  “ 

.  - 
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0  (o,o) 

0  (1,0)  0  (0,1) 

0  (2,0)  0  (1,1)  0  (0,2) 

0(3,0)  0(2,1)  0(1,2)  0(0,3) 

0  (4,0)  0  (3,1)  0  (2,2)  0  (1,3)  0  (0,4) 

0(5,o)  0(4,1)  0(3,2)  0(2,3)  0  (1,4)  0(0,5) 

0  (6,0)  0  (5,1)  0  (4,2)  0  (3,3)  0  (2,4)  0  (1,5)  0  (0,6) 


The  values  obtained  were 


1 


1  1 


73 


2  2  2 

5  6  6  6 

17  23  24  24  24 

109  118  120  120  120 

690  714  720  720 


388 


618 


720 
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Introducing 

(3.63)  u  =  ^(o,o)  +  /  (1,0)  j 

Cayley  used  (3*61)  and  (3*62) 

( 3  *  6U )  2  —  -  u  -  xu 

dx 

Hence  he  obtained 

g-x  +  lx2 

(3.65)  u  =  pzz.Tr— 

V  1  -  x 


♦  <f>  (2,0)  f 

to  show  that 

u 

1  -  x 


.  ..+  ^(m,0)  -  +  ...  , 

mi 


and 

(3.66)  <j>  ( m,0)  =  ml  f coefficient  of  x111  in 

Further  results  were 


(3.67) 


(m,l)  *  mi  (  coefficient  of  x111  in 

<^(m,2)  =  ml  |  coefficient  of  x01  in 

(f)(m,3)  -  ml  f  coefficient  of  x111  in 


hd^2/,,  3  ? 

e  /(I  -  x)  J 

2.  ei*  +  i*7(1  _ 

3J  eix  +  |xV(l-x 


From  (3.6I4)  5  Cayley  then  derived  the  following  relation 
for  finding  the  number  of  terms  in  a  symmetrical  determinant  of 
order  n 


(3.68)  (f>  (n,0)  =  n  <f>(n-l>0)  -  |(n-l)(n-2)  <f>  (n-3,0) 

We  shall  now  investigate  the  asymptotic  behavior  of  (m,n)  , 
for  large  m  .  Using  (3*61)  ,  we  find 


TOl-3 


i 


-  39 


(3.69)  2  jm  yi 


m,n=l 


mi  ni 


7 

y-i 


1  1  9 

2  X  +  l  X2 


1-x-y  V  1  -  x 


and 


(3.70) 


\ j>  (m,n) 

2  ri7---7  x111  y11 

mini 


m,n=o 


I— x 
1-x-y 


i 

2.  f 


V  1  -  X 


Hence 


(3.71)  ^  (m,n)  = 


dm 

c*>l 

1 _ 

(  1-x 

3  2  3C  +  k  X2  ) 

Sxm 

ro 

/  1-x-y  \f 

1-X  j 

y=0  . 

'  £ 

nl 

e  {  x  +  ix2 

dx™ 

(1  -  x)n+* 

x=0 

x=0 


and 


(3.72)  2  <Llh±_  jjn  =  n, 


m=o  ml 


i  l 


#  v2 


2  X  +  4X 


1  1  5 

IX  +  4  X2 


- - — ri  -  i  (-l)n  nl  - A 

(1  -  x)n  a  (x  -  l)n+z 


This  generating  function  has  a  branch  point  at  x=l  . 
the  notation  of  chapter  two. 


In 


(3.73)  g(x)  =  i  (-1)  nl  e 


i  ...  i 


4- 


and 


(3.7U)  g(r)(l)  =  i  (-Dn  nl  e 


|x2  +  x 

— .  e  * 

dot 


x-0 


If  we  make  the  substitution  z  =  —  i  x  3  t  =  -  i  . 

2  5 


(3.17)  yields 


. 


I 


-  ho  - 


(3.75) 


whence 


-r  (eix2  +  *) 
dxr 


H  ( 

r 


<- « (ii 


J  x=0 


(3.76)  g(r)(l) 


3  .  r 

i  (-l)n  nJ  e*  (!)  Hr(-  i)  . 


As  a  result  of  (3.21)  , 


[I]  2r"2s 


(3.77)  H  (-  i)  =  (-  i)1  r 1  2  , 

r  s=o  si  (r-2s) i 


md 


(3.78)  Ho(-  i) 
i) 

H2(-  i) 


-2i 


H3(-  i) 


20i 


Hence,  from  (2.33) 
(3.79)  {>  (m,n)  nj  nl  < 


■H) 


i  V  P  (m-r+n+  2  ) 


iP(n+  |  -r) 


H  (-  i) 


To  obtain  an  explicit  expression  for  <f>  (m,n)  ,  we  first 
apply  the  rule  of  Leibniz  to  (3.71)  •  i.e.. 


(3.80)  <f>  (m,n)  = 


nj 


=  nj 


r=c  W 

dr  ,  N-n« 

—  (1-x) 

cbF 

j 

"2 

' 

Jx=0L 

5  n 

r(n+r+  |  ) 

dm-r 

r=o  W 

F  (n+  l  ) 

<tem-r 

,m-r 


dxm“] 


i  i  ? 

(e5*  +  +X) 


1  1  p 

2  X  +  4  X2 


x=0 


0  • 


I 


, 


-  Ill  - 


1  . 


If  z  =  fix,  t  ■  -  |i  ,  (3.17)  becomes 

x11 


(3.81)  e*x+**2 

Hence,  using  (3.21)  , 

m— r  i  .  1  i 

(3.82)  1-  (e2*  * 

dx1"-1. 


"  / i\n  ,  i  ,  x“ 

2  2)  Hn(--> 

n=o  w  / 


1  \m~T  i 

H  C-f) 
m-r  2 


J  x=0 


(m-r) 


,  m 


>m-r 


s=o  si  (m-r-2s)l 


and 


(3.83)  </>(m,n)  = 


nl  ml  m  P(n+r+  \  )  2r  _ 

2  - . .  2 


ffi 


r=o  F'(n+ |  )  rl  s=o  si  (m-r-2s)i 


Let  us  now  derive  a  second  explicit  expression  for  (f)  (m,n) 
Substituting  x  =  u/n  in  (3.72)  ,  we  have 


(3.8U)  0(m,n)  - 


nl  n1 


du 


m 


(1 .  a  r 


u=0 


However, 


(3.85) 


u 

—  + 

i, 

(-)2 

u  +  i(u  + 

-2) 

*K{ 

’.a% 

■  a’) 

+  2-1 

n 

4 

In  1 

n 

2 

n2  \ 

2 

3  > 

r=3  nr  1 

,  2r 

r+l) 

(1  -  -  ) 
v  n  ' 


e"  1  t  2  A  (n )  •  - 
/  r=l  n 


where 


( 


; 


(3.86) 


-  U2  - 


A  (u) 


1  „ 

u  +  2  uz 


A2  (\x) 


A  (u) 


Au(u) 


= 

u2  + 

5  5 

6  U  + 

= 

5  ? 

I  u  + 

ll|u4 

7  5 

+  2+U 

1 

+  4*  U 

B 

6  5 
3U 

17  * 
3?u 

t 

+  16  ” 

It  can  be  shown  that  Ar(u)  is  a  polynomial  in  u  in 
which  the  least  power  of  u  is  r  ,  and  the  degree  is  2r  , 
Thus 


(3.37)  <f  (m,n)  = 


ni  n 


m 


,m  (  00 

—  e11  1+  2  A  (u)  - 

dura  l  r=l  r  nr 


u=0 


m 


=  nj  nw  1  +  Z 


r=l 


,m 


du 


eu  A  (u) 
m  r '  ' 


J  u=0 


Since 


d  '  sll  „S' 


(eu  u5) 


du 


ml 


J  u«0 


(m-s ) I 


j  we  obtain 


(3.88)  $  (m,n)  =  ni  nm  j  1  +  - 


+  +  JSL 

•  •  •  ~ 

m 


where 

(3.89) 


=  m  +  |  m(m-l) 

Bp  =  m(m-l)  +  §  m(m-l)  (m-2)  +  |  m(m-l)  (m-2)  (m-3) 

B  =  |  m(m-l)(m-2)  +  J-|  m(m-l)  (m-2)  (m-3)  +  ^.m(m-l) . . . (m-U) 


-^m(m-l).  ..(m-5) 


=  m(m-l)...(m-3)  +  §  m(m-l) . . .  (m-U)  +  i?m(m-l) . . .  (m-£) 

m 


**  ■  ‘  -  v  -■  3  '  Z6 
iLh  m(m-l)  * . .  (m-6)  +  3im(m-l). ..  (m-7) 


±  3 

38*  u 


1 


p 


,  * 


For  example, 


(3.90)  j  (3,n) 

<f)  (U,n) 

Example  Five  deals  with  the  number  f(n)  of  representations 
of  an  integer  n  as  an  ordered  product  of  factors  greater  than  1, 
The  asymptotic  behavior  of  f(n)  for  squarefree  n  was  determined 
by  Abe  Sklar  ,  He  defined  h(r)  =  f(Sr)  9  where  Sr 

represented  any  squarefree  integer  with  r  prime  factors,  and 
applied  the  Euler-MacLaurin  formula  to  the  result 

00 

(3.91)  h(r)  =  I  2  nr  2-n  , 

n=o 

to  obtain  the  following  theorem: 


-  nl  n3  /  1  ♦  i  ♦  il  ♦  I 


.  k  /  T  1°  17 

=  n  J  n  <  1  +  —  + 


2  n3 


7* 


(a)  For  every  positive  integral  r  , 

(3.92)  h(r)  =  2-1  rl  (log  2)~r~1  j  1  +  Rr  j  , 


and  the  remainder  R  can  be  exoressed  in  the  form: 
r 

(3.93)  2  2  (cos  @n)r+1  cos  |  (r-1)  6n  j  , 

n=l 


where  is  defined  by; 

n 


(3.9U)  cos  9 


log  3 
2rrn 


1  + 


e^r  j 


(b)  For  every  positive  integral  r  , 

r+1 


(3.9?) 


Rr  1  < 


^  Rr) 


where  f (r+1)  is  the  Riemann  zeta  function, 


■  ; 


- 


„ 

- 


i 


« 


■ 


' 


Let  us  first  develop  an  explicit  expression  for  h(r) 


Using  (3.91)  ,  we  find 


(3.96)  Z  h(r)  - 
r=o  r! 


\  2  2”n  enx 

n=o 


1  -  |ex 


Thus 


(3.97)  h(r)  = 


da* 


(1  - 


1  ex)_1 


x=0 


l 

=  2 


d* 

du 


(1  -  eu)’ 


u  =  -  log  2 


This  derivative  can  be  evaluated  using  the  Bernoulli  numbers, 
The  generating  function  of  these  numbers  is  given  by 


(3.98) 


eu  -  1 


GO  B 

„  mm 
Z  —  u 

m=o  ml 


whence 


(3.99) 


eu  -  1 


oo  B 


-  +  Z  — 

m=l  ml 


m  «-l 


Therefore, 


(3.100)  —  (eu  -  1)' 

dur 


(-l)r  ri 
„r+l 


um-r-X 


m=r+l  m  (m-r-1) .» 


(-DA  rl 

r+1 


”  Bk+r+l  uk 

2  — - - 

k=0  k+r+1  kj 


-  45  - 


and 


(3.101)  h(r)  = 


rl 


(log  2) 


r+1 


2 

k=o 


k+r+1 

k+r+1 


(-  log  ?-) 

kl 


k 


1 

2 


_j  !_Vl  (log  2)r+1 
r+1  )  r+i  rl 


(log  2) 

Br+2  (log  2)r+? 


+  ... 


r+2  r\  1J 

■where  the  Bernoulli  numbers  are  given  explicitly  by 


(3.102)  By 


_ 

n 

2 

m 

V 

Li 

(^)r  / 

\)  r11 

♦ 

m=l 

r=l 

m+1  V 

/ 

We  shall  now  derive  the  asymptotic  behavior  of  h(r)  from 
(3.101).  The  following  identity  will  be  useful. 


(3.103)  — 
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Let  us  first  suppose  that  r  =  2s  -  1  .  Then 


(3.10U)  Br+2  =  E^u  =  ...  =  B2n+1  -  o  , 


and,  using  (3.103)  , 

(3.105) 
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(3.108) 


Br+1  (log  2)5*1 
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Clearly,  all  the  terms  involving  even  Bernoulli  numbers 
will  be  asymptotic  to  zero. 

A  similar  result  follows  if  r  =  2s  . 

Hence,  we  have  shown 


(3.109)  h(r) 


r!  (log  2) 


-r-1 


The  above  result  can  be  obtained  very  easily  using  the 
method  of  Darboux.  From  (3.96)  ,  we  find  that  the  generating 
function  of  the  sequence  ^  h(r)  |  has  singularities  at 
x  =  log  2  +  2mri  .  Thus  the  nearest  singularity  is  at  x  =  log  2  . 
If  we  write  (3.96)  in  the  form 


(3.110)  2  h(r)  Si  =  - 1 -  x  ~  lpg  2 

r=o  ri  x  -  log  2  2  -  ex 

we  see  that  the  singularity  at  x  =  log  2  is  a  pole  of  order 
one  and  that,  in  the  notation  of  chapter  two, 

(3.111)  g(x)  =  (x  -  log  2)/(2  -  ex) 


In  this  case,  g(x)  has  a  removable  singularity  at  x  =  log  2  , 
and  can  be  made  regular  there  by  defining 

(3.112)  g(log  2)  =  lira  | 

x-log  2  2  -  ex 

Since  we  are  dealing  with  a  pole  of  order  one,  (2.33)  yields 
a  one-term  asymptotic  expansion  for  h(r)  ,  i.e., 


(3.113)  h(r)  ^  i  rl  (log  2) 
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The  study  of  the  asymptotic  behavior  of  sequences  occurs  in 
many  branches  of  mathematics.  The  solution  of  this  problem  by 
the  method  of  Darboux  requires  the  existence  of  a  generating 
function  with  a  finite  number  of  singularities  of  a  particular 
type.  Thus  a  glance  at  the  generating  function  determines 
whether  or  not  the  method  is  applicable.  Examples  three  and 
five  suggest  that  the  method  of  Darboux  is  not  as  well  known  as 
it  should  be. 

It  is  obvious  of  course ,  that  the  method  of  Darboux  does 
not  solve  the  most  general  case  of  the  asymptotic  behavior  of 
sequences.  For  example,  if  the  generating  function  is  an 
integral  function,  the  method  fails.  However,  several  other 
methods  have  been  developed  for  such  generating  functions.  One 
would  have  to  say  that  the  general  problem  of  the  asymptotic 
behavior  of  arbitrary  sequences  is  still  unsolved. 


. 

. 

- 


, 

r 


-  U8  - 


BIBLIOGRAPHY 


Poincare",  H.  Sur  les  integrales  irregulieres  des 
equations  lineaires.  Acta  Mathematica  8(1886)  , 
295  -  3hh  . 


Erdelyi,  A.  Asymptotic  Expansions.  Dover 
Publications,  Inc.  1956  . 

Darboux,  M.  G.  Memoire  sur  lf approximation  des 
fonctions  de  tres-grands  nombres,  et  sur  une  classe 
etendue  de  developpements  en  serie.  Journal  de 
Mathematiques  pures  et  appliquees  Series  3* 
li(l8?8),  5-56  and  377  -  Ul6  . 

Robinson,  R.  A  new  absolute  geometric  constant? 

The  American  Mathematical  Monthly  58(1951)  * 

)462  ~  U69  • 

Cayley,  A.  On  the  number  of  terms  in  a  symmetrical 
or  partially  symmetrical  determinant;  with  an  addition. 
Monthly  Notices  of  the  Royal  Astronomical  Society 
3)4(1873  ~  187U)  ,  303  -  307  and  335.  (The  Collected 
Mathematical  Papers  cf  Arthur  Cayley  9  ,  135  -  190  , 
Cambridge  ,  1896  •  ) 

Sklar,  A,  On  the  factorization  of  squarefree  integers. 
Proceedings  of  the  American  Mathematical  Society 
3(1952)  ,  701  -  705  . 


*  ■ 


- 


■ 


- 


r 


• .  <  ; 


■ 


* 


•  -  •  J  ’  •V>v¥*At  ■  ■■'x;4A:K 


